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Colloidal capsules can sustain an external osmotic pressure; however, for a sufficiently large pres- 
sure, they will ultimately buckle. This process can be strongly influenced by structural inhomo- 
geneities in the capsule shells. We explore how the time delay before the onset of buckling decreases 
as the shells are made more inhomogeneous; this behavior can be quantitatively understood by 
coupling shell theory with Darcy's law. In addition, we show that the shell inhomogeneity can 
dramatically change the folding pathway taken by a capsule after it buckles. 
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Many important natural or technological situations re- 
quire understanding thin, spherical shells; examples in- 
clude colloidal capsules for chemical encapsulation and 
release [1-3], biological cells [4, 5], pollen grains [6], sub- 
mersibles [7] , chemical storage tanks [8] , nuclear contain- 
ment shells [8], and even the earth's crust [9]. In many 
cases, the utility of such a shell critically depends on its 
response to an externally-imposed pressure. For small 
pressures, a homogeneous, spherical shell, characterized 
by a uniform thickness, supports a compressive stress, 
and it shrinks isotropically. Above a threshold pres- 
sure, however, this shrinkage becomes energetically pro- 
hibitive; instead, the shell buckles, reducing its volume 
by forming a localized indentation at a random position 
on its surface. For the case of a homogeneous shell, this 
threshold pressure can be calculated using a linearized 
analysis of shell theory [10, 11], while the exact morphol- 
ogy of the shell after it buckles requires a full nonlin- 
ear analysis [12-14]. However, many shells are inhomo- 
geneous, characterized by spatially-varying thicknesses 
and elastic constants [6, 15-18]. Such inhomogeneities 
can strongly influence how a shell buckles [6, 10, 19, 21- 
24] . Unfortunately, despite its common occurrence in real 
shells, exactly how inhomogeneity influences the onset of 
buckling, as well as the shell morphology after buckling, 
remains to be elucidated. A deeper understanding re- 
quires careful investigations of the buckling of spherical 
shells with tunable, well-defined, inhomogeneities. 

In this Letter, we use a combination of experiments, 
theory, and simulation to study the buckling of spheri- 
cal colloidal capsules with inhomogeneous shells of non- 
uniform thicknesses. We show that the onset of buckling, 
above a threshold external osmotic pressure, is well de- 
scribed by shell theory; however, even above this thresh- 
old, the capsules do not buckle immediately. We find that 
the time delay before the onset of buckling decreases as 
the shells are made more inhomogeneous; these dynamics 
can be quantitatively understood by coupling shell the- 
ory with Darcy's law for flow through a porous capsule 
shell, even for highly inhomogeneous shells. Moreover, 



we find that the shell inhomogeneity guides the folding 
pathway taken by a capsule during and after buckling. 
We use these insights to controllably create novel col- 
loidal structures using buckling. 

We fabricate monodisperse thin-shelled capsules us- 
ing water-in-oil-in- water (W/O/W) droplets prepared by 
micronuidics [25, 26]. The inner and outer phases are 
a 10 wt % solution of polyvinyl alcohol (PVA) of vis- 
cosity \i = 13.5 mPa-s, as measured using a strain- 
controlled rheometer, while the middle oil phase is a 
photo-polymerizable monomer solution. The PVA so- 
lution is less dense than the oil; as a result, after the 
droplets are produced and collected, the light inner wa- 
ter droplets gradually rise within them. This causes the 
oil to gradually thin on the top side of each droplet and 
thicken on the bottom [27]. We exploit this effect to pre- 
pare capsules, with outer radius Ro, and spatially- varying 
shell thickness h{0) ~ ho — ScosQ; is measured from the 
top of the gravitationally-oriented shell, ho is the average 
shell thickness, and 5 is the total distance moved by the 
inner droplet, as shown schematically in Fig. 1(a). The 
shell inhomogeneity can thus be quantified by the ratio 
5 /ho. We use UV light to polymerize the oil either as the 
capsules are produced in situ, or after different average 
waiting times, t w [28]; this enables us to prepare sepa- 
rate batches of capsules characterized by varying degrees 
of shell inhomogeneity [29, 30]. Some capsules are subse- 
quently washed in de- ionized water. The shell is a solid 
characterized by a Young's modulus E « 600 MPa [31]; 
importantly, while this shell is impermeable to Na + and 
CI - ions, it is permeable to water [32]. 

To probe their mechanical response, we subject in- 
homogeneous capsules, characterized by t w = 1 min, 
6/ ho ~ 0.2 and ho/Ro = 0.017, to an external os- 
motic pressure by injecting and gently mixing 20/iL of 
the capsule suspension into a fixed volume of NaCl so- 
lution, VnclCI ~ 130 — 400 /jL. We investigate the pres- 
sure dependence of buckling using NaCl concentrations 
in the range 0.063-2.165 M. Estimating the total volume 
of the injected capsules using optical microscopy allows 



us to calculate the final NaCl concentration of the outer 
phase, which then ranges from CNaCi — 0.055 — 2.068 M. 
These correspond to osmotic pressure differences across 
the shell of n = (2c NaC i + U out - U in ) x N A k B T = 
0.025 — 10.09 MPa, where Na is Avogadro's constant, 
fc# is Boltzmann's constant, T w 300 K, and IL ou t an d 
U in are the measured osmolarities of the fluids outside 
and inside the capsules, respectively, in the absence of 
NaCl. For each batch of capsules studied, we monitor an 
average of 75 capsules over time using optical microscopy. 

The osmotic pressure difference across these inhomoge- 
neous shells forces the capsules to buckle; we observe the 
abrupt formation of localized indentations in the shells, 
as shown in Fig. 1(b). For each osmotic pressure in- 
vestigated, the fraction of the capsules that buckle in- 
creases over time, eventually plateauing, as shown in Fig. 
1(c). We quantify this behavior by fitting this increase 
to an empirical exponential relationship, exemplified by 
the smooth lines in Fig. 1(c). The plateau value of this 
function yields a measure of the total fraction of the cap- 
sules that ultimately buckle over sufficiently long times, 
while the time constant of this function yields a measure 
of the time delay before the onset of buckling, r. For suf- 
ficiently large II, the total fraction of the capsules that 
ultimately buckle increases dramatically with increasing 
II, as shown by the grey circles in Fig. 1(d); this indicates 
that the capsules buckle above a threshold pressure, II* 
[33]. We empirically fit these data using the cumulative 
distribution function of a normal distribution, shown by 
the black line in Fig. 1(d); the mean value and standard 
deviation of this fit yield a measure of II* and the spread 
in II*, respectively. 

We study the geometry dependence of II* by per- 
forming additional measurements on inhomogeneous cap- 
sules, polymerized in situ, with different shell thick- 
nesses and radii; these are characterized by 5 /ho ~ 0.2, 
and ho/Ro = 0.019 or ho/Ro = 0.1. Similar to the 
ho/Ro = 0.017 case, for sufficiently large II, the total 
fraction of the capsules that ultimately buckle increases 
dramatically with increasing II, as shown by the red tri- 
angles and blue squares in Fig. 1(d). Interestingly, we 
find that the threshold buckling pressure II* ~ (ho/Ro) 2 
[Fig. 1(d), inset]; this observation is reminiscent of the 
prediction of shell theory for the buckling of a uniform 
shell [11], despite the fact that our capsules are inho- 
mogeneous. To understand this behavior, we consider 
the local deformability of an inhomogeneous shell at var- 
ious points on its surface. Because the 2D stretching and 
bending stiffnesses scale as ~ h and ~ h 3 [34], respec- 
tively, the thinnest part of the shell, where h ~ ho — S, 
should be the easiest to deform. We directly visualize 
that buckling begins at this "weak spot" using confocal 
microscopy of inhomogeneous capsules with fluorescent 
shells characterized by 6/ ho ~ 0.84 [Fig. 1(b), lower 
panel]. Consequently, we expect the onset of buckling 
to be governed by deformations in this part of the shell. 
To quantify this expectation, we apply shell theory to an 
inhomogeneous shell characterized by the same geome- 
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FIG. 1: (a) Schematic showing the capsule geometry investi- 
gated, (b) Upper: buckling of a capsule; scale bar is 20/xm. 
Lower: buckling begins at the thinnest part of the shell for 
capsules with thickness inhomogeneity 6/ ho « 0.84; scale bars 
are 50/xm. (c) Fraction of capsules buckled over time, for three 
different osmotic pressures II. Capsules have mean shell thick- 
ness ho = 1.2/im, outer radius Ro — 70/xm, and 6/ ho = 0.20. 
Smooth lines show exponential fits, (d) Total fraction of cap- 
sules that ultimately buckle over time for varying II, for cap- 
sules with ho,Ro, and 8 /ho — 1.2/im, 70/im, and 0.20 (grey 
circles), 1.3/im, 67/im, and 0.23 (red triangles), and 5.5/xm, 
55/im, and 0.19 (blue squares). Smooth curves are fits to the 
data using the cumulative distribution function of the normal 
distribution. Inset shows mean osmotic pressure of each fit 
versus ho/Ro] vertical and horizontal error bars show stan- 
dard deviation of each fit and estimated variation in ho/Ro, 
respectively. Straight line shows (ho/Ro) 2 scaling, (e) Time 
delay before the onset of buckling, r, normalized by for 
varying II, for the same capsules as in (d). Filled points show 
II > IT while open points show II < II*. Vertical error 
bars show uncertainty arising from estimated variation in ho- 
Black line shows II -1 scaling, (f) Time delay r decreases with 
the wait time before a shell is polymerized, t w ; capsules have 
ho = 1.2/im and Ro = 70/im, and are buckled at II « 0.86 
MPa > IT. Black line shows theoretical prediction coupling 
shell theory and Darcy's law, as described in the text, with 
k « 3.5 x 10" 24 m 2 . 



try as the experimental capsules [26] ; this analysis yields 

n * = 75ft 5 ) (^) 2 ~ 47 °(V#o) 2 MPa, assum- 
ing a Poisson ratio v ~ 1/3. The dependence of II* on 
ho~S confirms our expectation that the threshold buck- 
ling pressure is set by the thinnest part of the inhomoge- 
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neous shell. Moreover, we find U*(Ro/h ) 2 w 600 ± 200 
MPa for the experimental capsules [solid line, Fig. 1(d) 
inset] , in good agreement with our theoretical prediction. 
This indicates that the onset of capsule buckling is well 
described by shell theory. 

Within this framework, for II > II* , a capsule remains 
spherical before it buckles; it initially responds to the 
applied pressure by contracting uniformly, reducing its 
volume from its initial value, Vo, by a threshold amount 
A]/*, before buckling. We find that the time delay be- 
fore the onset of buckling, r, strongly decreases with in- 
creasing osmotic pressure II > II*, as shown by the filled 
points in Fig. 1(e). We hypothesize that this behavior 
reflects the dynamics of the fluid flow through the capsule 
shell; for the capsule to buckle, a volume AV* of fluid 
must be ejected from its interior. The time delay can 
then be estimated as r = AV*/Q [35], where both AV* 
and Q, the volumetric rate of fluid ejection from the cap- 
sule interior, are functions of 6/ ho. We calculate AV* for 
inhomogeneous shells using shell theory and validate the 
calculations with numerical simulations; the fluid ejec- 
tion rate Q follows from integrating Darcy's law over the 
surface of the capsule geometry shown in Fig. 1(a) [26]. 
Combining these results, we obtain 
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where Qo = 47i\RqIIA;/ fiho and k is the shell permeabil- 
ity. For the inhomogeneous capsules, characterized by 
6/ ho ~ 0.2, we thus expect t//iq ~ 0.8/i/HI; our experi- 
mental measurements of r allow a direct test of this pre- 
diction. Above II*, the data collapse when r is rescaled 
by ho, as shown by the filled points in Fig. 1(e), con- 
sistent with our expectation; moreover, by fitting these 
data [black line in Fig. 1(d)], we obtain an estimate for 
the shell permeability, k « 7 x 10 -24 m 2 . We use opti- 
cal microscopy to directly measure the rate at which the 
capsule volume decreases immediately after the onset of 
buckling [26]; this gives an independent measure of the 
shell permeability. We find k « 2 x 10" 24 m 2 [Fig. S8], 
in good agreement with the fit shown in Fig. 1(e); this 
further confirms the validity of Eq. 1. 

To test the applicability of this picture to even more 
inhomogeneous capsules, we measure r for capsules poly- 
merized at different t w ; these have shells with ho/Ro = 
0.017 and 6/ ho ranging from 0.2 up to 0.84. We impose a 
fixed osmotic pressure II « 0.86 MPa > II*. We observe 
that r decreases only slightly with increasing t w < 10 3 
s; however, as t w is increased above this value, r drops 
precipitously over one order of magnitude, as shown by 
the points in Fig. 1(f). To quantitatively compare these 
data to Eq. (1), we estimate the dependence of 6/ ho on 
t w using lubrication theory; we validate this calculation 
using direct measurements of 5 /ho for capsules prepared 
at varying t w [26, 36]. Remarkably, we find good agree- 
ment between our data and Eq. (1), with k w 3.5 x 10 -24 
m 2 , as shown by the black line in Fig. 1(f); in particular, 
this simple picture captures the strong decease in r at 
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FIG. 2: Folding pathways for different shell inhomogeneities. 
(a-c) Optical microscope images exemplifying buckling at 
II « 0.86 MPa of (a) slightly inhomogeneous capsules poly- 
merized in situ (t w « 0), with 6/ ho ~ 0.2, (b-c) very in- 
homogeneous capsules polymerized after a wait time t w — 1 
day, with 6 /ho ~ 0.84. Very inhomogeneous capsules buckle 
through the formation of either (b) one single indentation or 
(c) two indentations. At is time elapsed after buckling. Scale 
bars are 35/im. (d-e) Examples of simulated shells with sim- 
ilar geometries as the capsules shown in (a-c), for varying 
fractional volume reduction AV/Vo. Color scale indicates the 
spatially- varying shell thickness. 



t w ~ 10 3 s, with a shell permeability consistent with our 
independent measurements [Fig. S8]. While these results 
do not rule out other possible functional forms of r, they 
further suggest that the time delay before the onset of 
buckling can be understood by combining shell theory 
with Darcy's law for flow through the capsule shell, even 
for very inhomogeneous shells. 

The shell thickness inhomogeneity may continue to 
guide the development of deformations in a capsule af- 
ter it buckles. To explore this possibility, we use opti- 
cal microscopy to monitor the evolution of the capsule 
morphologies after the onset of buckling. Slightly inho- 
mogeneous capsules typically buckle through the sudden 
formation of a single circular indentation. As this inden- 
tation grows over time, its perimeter eventually sharpens 
into straight ridges connected by 2-3 vertices [14, 19, 20]; 
this folding pathway is exemplified by capsules polymer- 
ized in situ, characterized by 6/ ho ~ 0.2, as shown in Fig. 
2(a). This sharpening reflects the unique physics of thin 
shells: because it is more difficult to compress the capsule 
shell than it is to bend it, localizing compressive defor- 
mations only along sharp lines and points on the capsule 
surface requires less energy than uniformly compressing 
the shell [37]. Interesting differences arise for very inho- 
mogeneous capsules polymerized after t w = 1 day, char- 
acterized by 5 /ho ~ 0.84. The initial folding pathway 
is similar; however, the perimeters of the indentations 
formed in these capsules sharpen into straight ridges con- 
nected by 4-5 vertices, more than in the slightly inhomo- 
geneous case, as shown in Fig. 2(b). Moreover, surpris- 
ingly, roughly 30% of the very inhomogeneous capsules 
begin to buckle through the formation of one, then two, 
adjacent indentations, as exemplified in Fig. 2(c). The 
perimeters of these indentations grow over time, even- 
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FIG. 3: Colloidal capsules with two or three spherical in- 
terior compartments, schematized in left panels, buckle at 
"weak spots" (arrows). This forms shapes with two or three 
equally-spaced circular indentations after buckling (right pan- 
els). Scale bars are 100/im. 



tually meeting, coalescing, and sharpening into straight 
ridges connected by 4-5 vertices [Fig. 2(c)]. These obser- 
vations directly demonstrate that the deformations of a 
capsule after it buckles are sensitive to the shell inhomo- 
geneity. 

To gain insight into this behavior, we perform numeri- 
cal simulations of two different shells, a slightly inhomo- 
geneous shell with 5 /ho = 0.20, and a very inhomoge- 
neous shell with 5 /ho = 0.82, similar to the experimental 
capsules. As the shell volume is reduced below Vb — AV*, 
both shells buckle through the formation of a single in- 
dentation centered at the thinnest part of the shell, as 
shown in the leftmost panels of Fig. 2(d-e). As AV 
increases, this indentation grows and its edges sharpen. 
We find that the indentations formed in the very inhomo- 
geneous shells begin to sharpen at smaller AV/Vo, and 
ultimately develop more vertices than those formed in 
more homogeneous shells [Fig. 2(d-e)] [26]. These re- 
sults qualitatively agree with our experimental observa- 
tions [Fig. 2(a-c)], further confirming that after the onset 
of buckling, the folding pathway of a shell depends on the 
inhomogeneity. However, in contrast to the experimental 
capsules [Fig. 2(c)], we do not systematically observe the 
formation of adjacent indentations in the simulations on 
very inhomogeneous shells [38]. This presents a puzzle 
requiring further inquiry. 

Our capsules may be used to guide colloidal self- 



assembly; for example, a colloidal particle can sponta- 
neously bind to the indentation formed during buckling 
through a lock-and-key mechanism [39]. This mecha- 
nism is typically applied to a homogeneous colloidal 
particle, which buckles through the formation of a single 
indentation at a random position on its surface. We 
apply our findings to create multiply-indented capsules 
having two-fold or three-fold symmetry. To do this, we 
form double emulsions with two or three inner droplets 
of radii larger than half the radius of the outer droplet. 
Consequently, the inner droplets pack closely to form 
dimers or trimers [40], as shown schematically in Fig. 
3. The double emulsions are then polymerized, forming 
solid capsules with two or three spherical compartments 
in their interiors, and two or three equally-spaced "weak 
spots" in the capsule shell [arrows in Fig. 3]. When 
exposed to a sufficiently large osmotic pressure, these 
capsules buckle through the formation of multiple, 
equally-spaced indentations at the weak spots, as shown 
in Fig. 3. This approach is thus a versatile way to create 
capsules of desired symmetries, and extends the range 
of structures that can be used for lock-and-key colloidal 
assembly. 
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PACS numbers: 

I. MICROFLUIDIC FABRICATION OF 
CAPSULES 

To prepare monodisperse capsules, we use W/O/W 
double-emulsion droplets made with two different types 
of glass capillary microfluidic devices. To make cap- 
sules with ho/Ro = 0.017, we use a device composed 
of a hydrophobic tapered injection capillary (World Pre- 
cision Instruments, Inc., IB 100-6), pre-treated with n- 
octadecyltrimethoxyl silane (Aldrich), inserted in a sec- 
ond square capillary (AIT glass) [1]; the inner diame- 
ter of the square capillary and the outer diameter of 
the injection capillary are both 1 mm. Furthermore, 
a small tapered capillary is inserted into the injection 
capillary to simultaneously inject two immiscible flu- 
ids, as shown in Fig. SI (a). To confine the flow near 
the injection tip and thereby increase the flow veloc- 
ity, a hydrophilic circular capillary, pre-treated with 
2- [methoxy(polyethyleneoxy) propyl] trimethoxyl silane 
(Gelest, Inc.), is inserted into the square capillary at its 
other end. 

To make W/O/W double-emulsion droplets, we simul- 
taneously introduce two immiscible phases, 10 wt% aque- 
ous solution of poly (vinylalcohol) (PVA, Mw 13 000- 
23 000) and ethoxylated trimethylolpropane triacrylate 
(ETPTA) containing 0.2 wt% photoinitiator (2-hydroxy- 
2-methylpropiophenone, Aldrich), through the injection 
capillary at flow rates of 300 /iL/h and 150 /iL/h, respec- 
tively. The aqueous solution is saturated with ETPTA 
before injection to prevent diffusion of ETPTA molecules 
through water. We use the same aqueous solution of 
PVA as the continuous phase; this is injected through 
the square capillary at the flow rate of 3000 /iL/h. Be- 
cause of the hydrophobic nature of the injection capillary, 
the ETPTA flows along the inner surface of the injection 
capillary; by contrast, the aqueous solution flows through 
the center of the injection capillary as plug-like droplets 
that do not contact the wall. These plug-like droplets are 
emulsified at the tip of the injection capillary, resulting 
in monodisperse double-emulsion droplets with an ultra- 
thin middle layer as shown in Fig. Sl(b). The ETPTA 
formed between plug-like droplets produces large blobs at 
the tip of the injection capillary; we separate these from 
the double-emulsion droplets by exploiting their density 
difference. We prepare thin-shelled capsules by photo- 



polymerization of the droplet middle phase, ETPTA. 

To make capsules with ho/Ro — 0.019 and 0.1, we use 
a device composed of hydrophobic injection and collec- 
tion capillaries as shown in Fig. S2(a) [2]. The innermost 
aqueous phase is injected through the injection capillary 
and the middle ETPTA phase is injected through the in- 
terstices of the injection and the square capillaries. The 
continuous aqueous phase is injected through the inter- 
stices of the collection and the square capillaries as a 
counter flow to the innermost and the middle phases. 
These three streams flow coaxially through the orifice of 
the collection capillary, making double-emulsion drops in 
a dripping mode. For ho/Ro = 0.019, flow rates of the in- 
nermost, the middle, and the continuous phases are kept 
at 400, 80, 2000 /iL/h, respectively, as shown in Fig. 
S2(b). For ho/Ro — 0.1, flow rates of the innermost, the 
middle, and the continuous phases are kept at 200, 100, 
2000 /iL/h, respectively, as shown in Fig. S2(c). 

Capsules photopolymerized in situ are prepared in the 
following manner: As soon as the drops flow out from 
the nozzle of a microfluidic device to water bath contain- 
ing 8 wt% aqueous solution of PVA which is saturated 
with ETPTA, they are polymerized by continuous UV 
irradiation (Omnicure S1000). 

Capsules with t w = 5 s or 1 min are prepared in the fol- 
lowing manner: As they are continuously produced from 
the nozzle of the microfluidic device, the droplets are col- 
lected into a chamber containing 8 wt% aqueous solution 
of PVA saturated with ETPTA and are repetitively irra- 
diated with UV irradiation (Omnicure S1000) of duration 
1 second, fully polymerizing the ETPTA into a thin solid 
shell, every 2 x t w . This forms a population of capsules 
that have remained quiescent in the collection chamber 
for an average waiting time of t w before polymerization. 

Capsules with t w = 1 hr, 1 day, 2 days or 5 days are 
prepared in the following manner: As they are continu- 
ously produced from the nozzle of the microfluidic device, 
the droplets are collected for a total time of 10 min into 
a chamber containing 10 wt% aqueous solution of PVA 
saturated with ETPTA. They are then photopolymerized 
after a waiting time of t w — 5min. 

Optical microscope images of monodisperse microcap- 
sules with t w = 1 min are shown in Fig. S3. 

To make microcapsules with two or three spherical 
compartments, we employ two-step emulsification in a 
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capillary microfluidic device [3]. Inner water drops are 
produced at the tip of small tapered capillary in a drip- 
ping mode, which in turn are encapsulated into oil drops 
at the tip of injection capillary. Through precise control 
of flow rates of three streams, we can manipulate gener- 
ation frequency of inner water drops and outer oil drops. 
Therefore, the number of inner drops confined in outer 
drop is controllable. To make the thin shell, the flow rate 
of middle oil stream is maintained to be small enough to 
encapsulate two or three inner drops with nonspherical 
envelope. The nonspherical double-emulsion drops are 
polymerized by in-situ UV irradiation, resulting in mi- 
crocapsules with multiple spherical compartments. 

To buckle the capsules, we gently mix the capsule sus- 
pension with the NaCl solution by repeatedly aspiring 
and ejecting the fluid with a micro-pipette; this ensures 
the capsules are more directly exposed to the external 
osmotic pressure. 

We use either scanning electron microscopy (SEM) or 
confocal microscopy to directly measure the capsule ge- 
ometrical characteristics. 




FIG. SI: Microfluidic preparation of double-emulsion drops 
with ho/Ro = 0.017. Scale bar is 200/xm. 




FIG. S2: Microfluidic preparation of double-emulsion drops 
with (b) ho/Ro = 0.019 and (c) ho/Ro = 0.1. Scale bars are 
100/xm. 




FIG. S3: Monodisperse microcapsules. Scale bars in (a) and 
(b) are 500 jim. and 200 /mi, respectively. 



II. BUCKLING THEORY OF 
INHOMOGENEOUS SHELLS 



Before we derive the threshold buckling pressure II* 
and associated volume change AV* for inhomogeneous 
shells, it is instructive to study the derivation for homoge- 
neous shells of thickness ft, elastic modulus E, Poisson's 
ratio v and radius R. We use shallow-shell theory and fol- 
low the presentation by Hutchinson (ref. 4). The shallow- 
shell description involves isolating a shallow section of the 
shell of size ~ L and defining a Cartesian coordinate sys- 
tem (x, y) tangential to it (Fig. S5). "Shallowness" refers 
to the condition that the section is small compared to the 
radius, L/R <C 1, so that slopes of the surface are small; 
the validity of this approach for the buckling solution 
will be checked later. The shape of the middle surface 
is described by in-plane displacement fields v and a 
transverse displacement w of the middle surface of the 
shell. Ignoring terms of order (x/R) 2 , (y/R) 2 <C 1, the 
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FIG. S4: SEM showing buckled shell with inhomogeneity. 
Shells buckle at the thinnest part; as a result, the thinnest 
part contacts the thickest part of the shell. Right panel: the 
left arrow shows the thickest part of the shell, and the right 
arrow shows the thinnest part of the shell; the solid in be- 
tween the two parts of the shell is precipitated PVA. Scale 
bar is 2/xm. 




FIG. S5: Coordinates and displacements in shallow shell the- 
ory. 



displacement-strain relations for the section are 
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The elastic energy of the shell involves stretching, 
bending and a pressure term: 
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where the stresses are related to the strains via 
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The initial response of the shell to an external pressure 
II is to contract uniformly, building up uniform compres- 
sive shell stresses 



'- J xx l - J yy TLR/2 

with corresponding transverse displacement field 
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that is also uniform over the whole section. The elas- 
tic energy associated with this compressed state, which 
does not involve bending of the shell surface, is roughly 
a 2 /(Eh) - U 2 R 2 /(Eh) per unit area. When II gets 
large, it becomes energetically favorable for the shell 
to introduce some bending deformations to trade the 
elastic energy of pure compression for bending energy. 
The buckling transition is identified by the value of II 
for which a nonuniform deformation mode of the shell 
first arises which is energetically favorable compared to 
the uniform compression. A linear buckling analysis of 
the Euler-Lagrange equations associated with the en- 
ergy functional (Eq. S2) shows that a nontrivial solution 
for the transverse displacement field w first arises at a 
threshold buckling pressure [4] 
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The form of the associated displacement field is si- 
nusoidally, w — Ae zqx , with period set by a two- 
dimensional wavevector q that satisfies |q| = (f = 



[12(1 



(The corresponding solutions for 



the in-plane displacement fields u and v have similar 
forms.) Since a multitude of two-dimensional wavevec- 
tors exists with magnitude equal to g*, there are many 
degenerate buckling modes at the transition. Crucially, 
the wavelength associated with the buckling modes is ex- 
tremely small for thin shells: 



A* = 27r/g* 
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This observation justifies the use of shallow-shell theory 
to calculate the buckling pressure — a shallow section of 
the shell with L/R <C 1 can still accommodate many 
wavelengths of the buckling modes and provides a good 
description of the elastic energy of these modes. 

We now consider an inhomogeneous shell with a con- 
tinually varying shell thickness h(6) = ho — 5 cos 6 and 
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radius R of the middle surface. By analogy with the 
buckling of homogeneous shells, we expect the buckling 
transition to be signalled by the existence of a nontriv- 
ial solution to the energy-minimizing displacement fields 
of the middle surface. Furthermore, we expect such a 
mode to first appear in the vicinity of the thinnest point 
of the sphere, = 0. We consider a shallow section 
of the sphere, with Cartesian coordinates whose origin 
coincides with the thinnest point. As before, shallow- 
ness implies that the section is small enough that slopes 
are small compared to the base of the section. The 
elastic energy and the stress-strain relations are modi- 
fied because the thickness of the shell is now position- 
dependent: h(x,y) = ho — S (1 — Z(x, y)/R) , where 
Z(x, y) = \/R 2 — x 2 — y 2 is the shape of the undeformed 
middle surface of the shell. Expanding in powers of x/R 
and y/R, we have 



h(x,y) = h - S 
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but the shallowness of the section implies that 
(x/R) 2 , (y/R) 2 <C 1 and the section is well-described by 
setting h « ho — 5, a constant. Thus, the only modifica- 
tion to the shallow-shell theory calculation of the critical 
buckling pressure of the uniform shell is to set the shell 
thickness toh = ho — S, from which we readily obtain the 
buckling pressure 
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The wavelength of the buckling modes that arise at this 
pressure is A* = 2tt v / (/i - 6) R/ [12(1 - v 2 )} 1 ^ <C R, 
justifying the assumptions made. Essentially, the spatial 
thickness variation of the capsules under study is so gen- 
tle that there is hardly any variation in the mechanical 
properties of the capsule at the scale of the small wave- 
length associated with the buckling modes. Thus, the 
buckling pressure is set by the thickness of the capsule 
in the vicinity of the thinnest spot, where it is approxi- 
mately ho — 5. (We could consider different sections of 
the capsule, with higher thicknesses, but the associated 
buckling mode would arise at a pressure higher than of 
Eq. S9.) Considering the magnitude of the terms ex- 
cluded from the Euler-Lagrange equations by our approx- 
imation of constant thickness shows that the corrections 
to Eq. S9 from taking into account the spatial variation 
in thickness over the shallow shell would be smaller by a 
factor 0(5/R). 

To test our theoretical arguments, we perform numer- 
ical simulations on shells with 0.01 < ho/R < 0.04 and 
< 5 /ho < 0.9 (see section VI below for details). The 
results are reported in Fig. S6(a), showing good agree- 
ment between the theoretical prediction and the simula- 
tion results. The slight discrepancy between theory and 
simulation likely reflects the sensitivity of the buckling 
transition to non-uniformities in the simulation mesh. 
This sensitivity is due to nonlinear couplings between 



the degenerate buckling modes that arise as the classi- 
cally predicted buckling pressure is approached, which 
drive the transition to happen below the buckling pres- 
sure predicted by the linear stability analysis [4]. We 
expect this sensitivity to be more pronounced for the 
uniform case, where the unstable modes extend over the 
entire shell, in contrast to the nonuniform case, where 
the unstable modes are localized in the shallow region 
near the thinnest part of the shell. 

We now turn to the threshold volume change at buck- 
ling, Ay* . The volume change prior to buckling is deter- 
mined by the initial prebuckling response of the shell to 
the external pressure. For a spherical shell, the only ax- 
isymmetric stress distribution in response to a uniform 
pressure p that is well-behaved at the poles is the uni- 
form stress distribution, Eq. (S4) [5]. In contrast to the 
uniform shell, however, the resulting transverse displace- 
ment varies with the shell thickness: 
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The volume change in response to pressures up to the 
buckling pressure is thus (to lowest order in the inward 
displacement) 
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The threshold volume change immediately before buck- 
ling is obtained by using II* (Eqn. S9) in the above ex- 
pression to get 
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where Vo = 47ri? 3 /3 is the initial volume of the shell. 
We recover the result for uniform shells, AV* /Vo = 
y/Z(l-v)/(l + v) x ho/R (Ref. 24) in the limit 5 /ho 
0. Fig. S6(b) compares the analytical expression to nu- 
merical simulations on shells with various thicknesses and 
inhomogeneities, showing good agreement. Moreover, we 
find that the total fraction of capsules that ultimately 
buckle in our experiments increases with increasing t w , 
and hence increasing inhomogeneity 5 /ho, as shown in 
Fig. S7. This supports the theoretical finding that inho- 
mogeneity reduces the mechanical strength of the shell. 

We note that the linearized buckling analysis sketched 
out here can describe the shell up to the point of the 
buckling transition, but cannot describe the postbuck- 
ling shape of the shell. Once the buckling transition is 
reached, the nontrivial high- wavelength buckling mode, 
when it arises, coalesces into a single inversion (as shown, 
for example, in [10]) unless an additional mechanism is 
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introduced to arrest this coalescence for example, by 
placing a slightly smaller, rigid, concentric sphere within 
the shell [26]. In Carlson et aVs work, the inner sphere 
prevented the unstable buckling mode from coalescing 
into a single inversion, giving rise instead to a regular 
dimple pattern. In the absence of such a stabilizing mech- 
anism, the high- wavelength mode itself is not observed in 
buckling experiments including the ones on the inhomo- 
geneous capsules prepared here. 



III. FLOW RATE OUT OF AN 
INHOMOGENEOUS CAPSULE 



The buckling of a capsule is driven by the fluid ejec- 
tion from the capsule interior, due to the imposed osmotic 
pressure difference across the capsule shell, II. This is re- 
sisted by the mechanical pressure required to compress 
the shell, U rn , at most U m « 2E(1 - - 
1/Ro) w 3E(ho/Ro)[(l - AV*/V )~ 1/3 - 1], where R* is 
the radius of the shell at buckling. For the shells studied 
in this work, II m < II*; we thus expect the buckling dy- 
namics to be dominated by the imposed osmotic pressure 
for II ^> II* , and we neglect II m in the simple model pre- 
sented here. A full treatment of the buckling dynamics 
that explicitly includes II m confirms the validity of this 
simplification for the range of II studied here [7] . 

We now estimate the flow rate out of an inhomogeneous 
capsule due to II; the capsule has shell thickness h{6) = 
ho — 5cos9, as shown in Fig. 1(a) of the main text. We 
use spherical coordinates (r, 0, <p) centered on the capsule 
center. For an arbitrary area element dA on the shell 
surface at (i?o, 0, (/>), the local volumetric ejection rate is 
given by Darcy's law, dA ■ TLk/ jah(0)^ where k is the shell 
permeability and \i is the fluid viscosity. Integrating this 
over the entire shell surface yields the total ejection rate 
through the shell: 
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The time delay before the onset of buckling, r, is the 
time taken for the volume of fluid ejected from the shell 
to equate to the threshold buckling volume, r = AV* /Q. 
From Eq. S12 and Eq. S13, we obtain 
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FIG. S6: (a) Effect of inhomogeneity on the buckling pres- 
sure. The symbols show the buckling pressure from simu- 
lations for shells with various average thickness ho and in- 
homogeneity 8, normalized by the corresponding buckling 
pressure for uniform shells with the same average thickness, 
IIJUo = 2E/y/3(l - v 2 ) x (h /R) 2 (Eq. S6). The blue line 
shows the theoretical prediction, Eq. S9, and the yellow line 
shows 90% of the theoretical value, (b) Threshold volume 
change at buckling, from simulations. The blue line shows the 
theoretical prediction (Eq. S12), with v — 1/3, while the yel- 
low line shows 90% of the theoretical value. In both (a) and 
(b), the theory agrees with the simulations to within about 
90%, and the theory effectively captures the dependence of 
the mechanical strength on the inhomogeneity. The system- 
atic discrepancy of 10% or so (higher for the extremely thin 
shells with ho/R = 0.01) is likely a result of the sensitivity 
of the buckling transition to the small amount of disorder in 
the underlying mesh. It is known that the buckling transition 
of spherical shells is highly sensitive to imperfections in the 
shell [4] which significantly reduce the buckling pressure and 
through it the associated volume change. 



^ V / 3(1 - v) h ( S \ 2 

where Vb = 47r.Ro/ 3 ^ s ^ ne initial capsule volume, Qo = 
AirRyik/ iiho and k is the shell permeability. This is Eq. 
(2) in the main text. 



IV. ESTIMATE OF CAPSULE PERMEABILITY 

We estimate the capsule permeability k by measur- 
ing the change in the radius of a circular indentation r 
over time, immediately after it is formed in the shell, 
as shown in Fig. S8. We assume the volume of the in- 
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FIG. S7: Total fraction of capsules that ultimately buckle 
increases with wait time before polymerization t w , and hence 
the shell inhomogeneity. Capsules have ho = 1.2/xm and Ro = 
70/im, and are buckled at II « 0.86 MPa> II*. 



dentation AV cap is given by that of a spherical cap hav- 
ing radius r(t). The permeability is then given by k ~ 
^h d(AVca P yd(At) ^ reality, the edges of the indentation 

are rounded, with radius of curvature ~ ^/hoRo ~ 15/im 
[8]; we thus expect our estimated AV cap to overpredict 
the actual AV and consequently, we expect to underpre- 
dict the shell permeability by a factor ~ 2. 



V. LUBRICATION THEORY OF THE 
FORMATION OF THE INHOMOGENEOUS 
SHELL (COMMUNICATED BY HOWARD A. 
STONE, FEBRUARY 2012) 

To estimate the variation of the shell inhomogeneity, 
5 /ho, with time, t w , we apply lubrication theory to the 
double emulsion geometry shown in Fig. 1(a) of the main 
text. Before UV polymerization at time t w , the inner 
droplet containing 10wt% PVA has radius Ro — ho and 
the outer droplet containing ETPTA monomer has radius 
Ro] the shell thickness is then h(t w ) = ho — 5(t w )cos9. 
The droplets are collected in 8 wt% or 10 wt% PVA so- 
lution. We note that, for our shells, ho/Ro < 0.1. More- 
over, using our experimental measurements of S(t w ), we 
estimate the characteristic translation speed of the in- 
ner droplet as ~ 1/im/day; using a characteristic droplet 
length scale < 100/im, shell viscosity \i Q = 65mPa.s, 
fluid- fluid surface tension ~ 2.5mN/m, this corresponds 
to a Reynolds number Re < 10 -10 and a Capillary num- 
ber Ca < 10 -10 . The fluid- fluid interface can thus be 
approximated as spherical, and the flow in the shell be- 
tween the inner and outer droplets is well described using 
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FIG. S8: Change in the volume of a circular indentation 
formed in a capsule, AV cap , over time At. The indentation 
forms at At = 0. Top left panel shows the top view of the in- 
dentation formation, measured using optical microscopy; we 
use image processing to detect the edge of the indentation, 
shown in the panel to the right, and track the radius of the 
indentation over time, r(t). We assume a spherical cap ge- 
ometry, schematically shown in the side view, to calculate 
AVcap ; the data are shown for three different capsules (dif- 
ferent colors). We fit the small-time dynamics (At < 10 2 s) 
to measure the permeability. The capsules have ho/Ro = 0.1 
and 5 /ho « 0.2, and are buckled at II = 10 MPa. 



lubrication theory. We denote x — rO as the direction 
along the shell, where r is the spherical radial distance 
measured from the center of the inner sphere, and z is 
across the shell. 

Because both inner and outer droplets are stabilized by 
8-10 wt% PVA, Marangoni stresses resist shear stresses 
at the fluid interfaces; thus, we assume no-slip boundary 
conditions on the fluid-fluid interfaces at z = and z = 
h(6,t w ). The Stokes equation then yields the velocity 
distribution: 



dp 



2fi(Ro - ho) 80 



z(z-h(0,t w )) (S15) 



Integrating the continuity equation twice across the 
gap, with u r = at the outer boundary, u r — UcosQ at 
the inner translating boundary, and uq = u x = at both 
boundaries, we find the pressure distribution 



p(0,t w ) = po(t w ) 
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Balancing forces in the z direction, 
27r/Jl pn-e a sin<9d<9 = F\ where F b = 
(4/3)7r(i?o — ho) 3 Apg is the buoyant force on the 
inner sphere and Ap is the difference in density between 
the inner and outer spheres (0.08 g/cm 3 ). We complete 
the integral and substitute U = dS/dt W: and take the 
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limit 6/ ho — > 1 [25]; we obtain 
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To test this prediction, we use SEM or confocal mi- 
croscopy measurements of capsules with fluorescently- 
labeled shells to directly measure 8 /ho for capsules of 
varying t w . We find good agreement between the two, 
particularly for large t w , when 5 /ho precipitously rises to 
1, as shown in Fig. S9. We find 6/ ho — » 0.2 as t w — >> 0, 
in contrast to the theory; this is to be expected, due to 
the slight bias in the inner droplet position as the double 
emulsions are formed in the microfluidics device. 
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FIG. S9: Measurements of 5 /ho as a function of t w - Qfi ( R _ ho ^ , 
along with lubrication theory solution (black line), showing 
agreement between the two. 
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VI. NUMERICAL SIMULATIONS OF THIN 
INHOMOGENEOUS SHELLS 

To test our theoretical predictions and to understand 
the postbuckling shapes of the capsules, we perform nu- 
merical simulations of a variety of thin, elastic, spheri- 
cal shells whose surfaces are composed of "amorphous" 
randomly-positioned points connected by springs. Simi- 
lar models have been used to study the deformations of 



uniform spherical shells [9, 10], viruses [11-15] and pollen 
grains [16]. 

We simulate shells with a non-uniform thickness pro- 
file, described by h(6) = ho — ScosO, similar to the ex- 
perimental capsules; is the inclination angle from the 
thinnest point of the shell. The spring stiffnesses are 
chosen so that the 3D Young's modulus of the shell is 
E = 600 MPa and the Poisson ratio is v — 1/3, sim- 
ilar to the material making up the polymeric capsules. 
Since h is always much smaller than the radius of the 
shell, the elastic energy may be described by a sum of 
bending and stretching energies of an initially spherical 
two-dimensional elastic sheet [8, 17]: 
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where is the 2x2 strain tensor, kij is the change in the 
curvature tensor from its initial value, and the integration 
is carried out over the mid-surface of the shell. The initial 
curvature tensor k% of the shell mid-surface is that of a 
sphere with radius r: k^ = k\ 2 = 1/r; fcf 2 = at every 
point for an orthonormal basis set up in the tangent plane 
to the sphere at that point. 

To numerically simulate the shells, we discretize the 
elastic energies on a mesh of 20,000 points. The ini- 
tial, unstrained configuration is obtained by distributing 
the points quasi-randomly on the surface of a sphere, 
maintaining a minimum distance between nearest neigh- 
bors. The disorder in the initial mesh eliminates the 
effect of the 12 regularly spaced five- fold disclinations 
that inevitably arise when attempting to cover a spher- 
ical surface with equilateral triangles [11]. Bonds are 
drawn between nearest neighbor pairs to form a mesh; 
the topology of the mesh is unchanged after initializa- 
tion (no rearrangement of bonds takes place) , consistent 
with the polymerization imposed in the experiments. A 
point is arbitrarily chosen to be the thinnest point of 
the sphere. The elastic stretching energy of deformations 
from the initial unstrained configuration is approximated 
by a harmonic spring energy associated with each bond 
[18]: 
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(ij) 



where and rj- are the lengths in the deformed and ini- 
tial states of the bond connecting nearest-neighbor mesh 
points i and j, and Oij is the inclination angle between the 
centre of the spring and the thinnest point. The spring 
constant is chosen to reproduce the thickness-dependent 
elastic modulus in the continuum limit, with v = 1/3 
[18]. 
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The often-used discretization of the bending energy in 
terms of the angles between normals of adjacent facets 
in the mesh is not suitable for disordered meshes where 
all facets are not equilateral triangles [19]. Furthermore, 
this discretization scheme does not allow for a nonzero 
background curvature. Here we use a different approach. 
The curvature tensor is constructed from the mean cur- 
vature H and Gaussian curvature K (respectively, the 
trace and the determinant of the curvature tensor) as- 
sociated with each point in the deformed state. Note 
that the Gaussian curvature energy is not an invariant 
for inhomogeneous shells, and thus cannot be ignored, in 
contrast to uniform shells [20] . An approximation for the 
mean curvature associated with site i is [19]: 
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where <j^- = 



; (cot#i + cot#2)/2 is the length of a bond 
in the dual mesh connecting the centers of the two tri- 
angular facets ijk and ijk' that share the bond linking 
points i and j (calculated in terms of the interior angles 
61,62 of the vertices k and fc'), &i = Ylj(i) a ij r ij/^ is the 
area of the dual cell (the polygon of bonds of the dual 
mesh surrounding vertex z), is the three-dimensional 
coordinate vector of point z, and all sums are over the 
neighbors of site i. The Gaussian curvature is approxi- 
mated by exploiting the Gauss-Bonnet theorem [21]: 
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where aj is the angle facet j subtends at point i and the 
sum runs over all facets sharing point i. The discretized 
version of the bending contribution to the elastic energy 
[Eq. S3] is 
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We add a penalizing term E vo \ = A(V — Vr) 2 to the en- 
ergy, where V is the volume enclosed by all the facets of 
the mesh, with A initially chosen to be very small. We 
increment A systematically between successive minimiza- 
tions of the total energy, until the difference between the 
actual and target volumes is negligible. The final config- 
uration is checked for stability to a small random pertur- 
bation of all the point positions. 

To find the threshold volume reduction at collapse, we 
initially set Vp to the volume Vo enclosed by a sphere 
of radius r. We then reduce the target volume by small 
amounts, finding the minimum energy configuration after 
each decrement. The volume decrement between steps 
is 0. 00025 Vo. To quantify the buckling, we define the 
asphericity of the shell as [11] 
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where Ri is the radial distance of point i and (...) de- 
notes averaging over all points of the mesh. The discrete 
nature of the buckling inversion can be captured by plot- 
ting the asphericity of the final shape against AV/Vb- 
the asphericity first becomes nonzero at a finite volume 
reduction which we identify as AV* . 



A. Postbuckling shapes of shells 



The total energy E s 



■E b ,a, function of the 3 x 20, 000 
variables describing the point positions, is numerically 
minimized using the BFGS quasi-Newton optimization 
algorithm implemented in the freely available GNU Sci- 
entific Library [22]. 

To numerically obtain the threshold buckling pressure 
from simulations, we add a term E p = HV to the dis- 
cretized elastic energy energy, where V is the volume 
enclosed by all the facets of the mesh. The pressure II 
is gradually increased until the shell collapses, and the 
corresponding threshold pressure II* identified. 

To find energy-minimizing configurations at a particu- 
lar target volume Vr, we use the penalty method [16, 23]. 



To investigate the influence of inhomogeneity on the 
shape of the inversion at large volume reductions, we sim- 
ulate large volume reductions for two sets of geometric 
parameters: shells with ho/R = 0.018 and 6/ ho = 0.20, 
similar to the more homogeneous shells described in the 
main text, and shells with ho/R = 0.0165 and 6/ ho = 
0.82, similar to the very inhomogeneous shells. We use 
the penalty method outlined in the previous subsection 
to reduce the enclosed volume in steps of 0.0025Vb from 
Vr = Vo to Vr = 0.75Vb- For each set of geometric pa- 
rameters, we investigate the robustness of the resulting 
shapes by performing eight different simulations, in each 
case varying the point on the underlying mesh that corre- 
sponds to the thinnest point of the shell. This effectively 
(S22)changes the random discretization of the shell in each 
instance. 

In all instances, the shells buckle when the shell vol- 
ume is reduced past AV* (whose value is approximately 
0.01 15 Vo for the more homogeneous shells and 0.0008Vb 
for the very inhomogeneous shells) to form a single circu- 
lar indentation centered at the thinnest part of the shell. 
At larger volume changes, however, the shape of the in- 
version differs between the two shells and among different 
simulation instances for the same shell. In more homoge- 
neous shells, the inversion remains circular upto a volume 
change of AV ~ 0.05Vb, beyond which the perimeter 
sharpens into straight ridges connected by 4-5 vertices 
[main text, Fig. 3(d)]. Three out of the eight simulation 
runs on the more homogeneous shells displayed inversions 
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with 4 vertices, while the remaining five runs displayed 
inversions with 5 vertices. Remarkably, the two types 
of inversions do not differ significantly in their elastic 
energies, which explains why slight differences in the un- 
derlying mesh have a significant effect on the final shape. 

The shape evolution of the inversions in very inhomo- 
geneous shells is qualitatively different [main text, Fig. 
3(e)]. The inversion takes on a ridged appearance almost 
immediately after its formation, at a much smaller vol- 
ume change (AV ~ 0.004Vb) compared to the more ho- 
mogeneous shells. The number of vertices in the perime- 
ter is also consistently higher for the very inhomogeneous 
shells: out of eight simulation runs, three displayed in- 
versions with 5 vertices, four had inversions with 6 ver- 
tices, and one run displayed an inversion with 7 vertices. 
Again, there is no significant difference in the elastic en- 
ergies of the competing shapes. 

We note that while the exact number of vertices formed 
differs between the experiments and the numerical sim- 
ulations, both show that the inversions in very inhomo- 
geneous capsules develop more vertices than the more 
homogeneous case. 

We further characterize the dependence of postbuck- 
ling shapes on shell inhomogeneity by simulating large 
volume reductions for shells with ho/R = 0.2, similar to 
the experimental capsules, but varying 6/ ho from 0.05 to 
0.95. The number of vertices observed in the inversion is 
shown in Fig. S10. We find that the number of vertices 
increases with inhomogeneity and with increasing volume 
reduction. 



VII. MICROGRAPHS OF BUCKLED 
CAPSULES 
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scale bar is 20/im. 
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FIG. S12: Buckling begins at the thinnest part of the shell 
for capsules with thickness inhomogeneity 6/ ho ~ 0.84 
; scale bars are 50/im. These capsules are prepared with 
fluorescent shells to enable visualization with confocal 
microscopy. 




FIG. S10: Inhomogeneity dependence of the postbuckling 
shape. The number of vertices, n, of the polygonal inver- 
sion is shown for four volume reductions as a function of the 
shell inhomogeneity 6/ ho. A value of n — corresponds to 
a roughly circular inversion with no discernible sharp points. 
For all shells, ho/R = 0.02. The volume was reduced in steps 
of 0.005V& from V T = Vo to V T = O.S7V . We used the same 
mesh for all the simulations. 
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FIG. S13: Folding pathways for different shell inhomogeneities. (a-c) Optical microscope images exemplifying buckling at 
II « 0.86 MPa of (a) slightly inhomogeneous capsules polymerized in situ (t w ~ 0), with 6/ ho ~ 0.2, (b-c) very inhomogeneous 
capsules polymerized after a wait time t w = 1 day, with 6/ ho « 0.84. Very inhomogeneous capsules buckle through the 
formation of either (b) one single indentation or (c) two indentations. At is time elapsed after buckling. Scale bars are 35/im. 
(d-e) Examples of simulated shells with similar geometries as the capsules shown in (a-c), for varying fractional volume reduction 
AV/Vq. Color scale indicates the spatially- varying shell thickness. 



